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ABSTRACT 

Using numerical methods, we systematically study in the framework of ideal MHD the effect of mag- 
netic fields on heat transfer within a turbulent gas. We measure the rates of passive scalar diffusion 
within magnetized fluids and make the comparisons a) between MHD and hydro simulations, b) between 
different MHD runs with different values of the external magnetic field (up to the energy equipartition 
value), c) between thermal conductivities parallel and perpendicular to magnetic field. We do not find 
apparent suppression of diffusion rates by the presence of magnetic fields, which implies that magnetic 
fields do not suppress heat diffusion by turbulent motions. 

Subject headings: turbulence - ISM: general - galaxies: clusters: general - MHD 



1. ASTROPHYSICAL MOTIVATION 

It is well known that Astrophysical fluids are turbulent 
and that magnetic fields are dynamically important. One 
characteristic of the medium that magnetic fields and tur- 
bulence may substantially change is the heat transfer. 

There are many instances when heat transfer through 
thermal conductivity is important. For instance, thermal 
conductivity is essential in rarefied gases where radiative 
heat transfer is suppressed. This is exactly the situation 
that is present in clusters of galaxies. It is widely accepted 
that ubiquitous X-ray emission due to hot gas in clusters of 
galaxies should cool significant amounts of the intraclus- 
ter medium (ICM) and this must result in cooling flows 
(Fabian 1994). However, observations do not support the 
evidence for the cool gas (see Fabian et al. 2001) which 
is suggestive of the existence of heating that replenishes 
the energy lost via X-ray emission. Heat transfer from the 
outer hot regions can do the job, provided that the heat 
transfer is sufficiently efficient. 

Gas in clusters of galaxies is magnetized and the conven- 
tional wisdom suggests that the magnetic fields strongly 
suppress thermal conduction perpendicular to their direc- 
tion. Realistic magnetic fields are turbulent and the issue 
of the thermal conduction in such a situation has been 
long debated. A recent paper by Narayan & Medvedev 
(2001) obtained estimates for the thermal conductivity of 
turbulent magnetic fields, but those estimates happen to 
be too low to explain the absence of cooling fiows for many 
of the clusters of galaxies (Zakamska & Narayan 2002). 

Narayan & Medvedev (2001) treat the turbulent mag- 
netic fields as static. In hydrodynamical turbulence it is 
possible to neglect plasma turbulent motions only when 
the diffusion of electrons which is the product of the elec- 
tron thermal velocity Veiect and the electron mean free 
path in plasma Imfp, i-e. Veiectlmfp, is greater than the 
turbulent velocity Vturb times the turbulent injection scale 
linjj i.e. vturbknj- If such Scaling estimates are applica- 



ble to heat transport in magnetized plasma, the turbu- 
lent heat transport should be accounted for heat trans- 
fer within clusters of galaxies. Indeed, data for Veiectlmfp 
given in Zakamska & Narayan (2002; Narayan & Medvedev 
2001) provide the classical Spitzer (1962) diffusion coeffi- 
cient Ksp = Veiectlmfp ~ 6.2 X 10^*^ cm^ scc"^ for the inner 
region of i? ~ lOOkpc and Ksp = Veiectlmfp ^ 3.6 x 10^^ 
cm^ sec^^ for the very inner region of ~ lOkpc (for 
Hydra A). If turbulence in the cluster of galaxies is of 
the order of the velocity dispersion of galaxies, while the 
injection scale is of the order of 20 kpc, the diffusion co- 
efficient is ~ Vturbknj ^ 3.1 X lO'^" cm^ sec~^, where we 
take Vturb 500 km/sec. 

Earlier numerical studies by Cho, Lazarian & Vishniac 
(2002) revealed a good correspondence between hydrody- 
namic motions and motions of fluid perpendicular to the 
local direction of magnetic field. To what extend heat 
transfer in a turbulent medium is affected by a magnetic 
field is the subject of the present study. To solve this prob- 
lem we shall systematically study the passive scalar diffu- 
sion in a magnetized turbulent medium, compare results 
of MHD and hydrodynamic calculations, and investigate 
the heat transfer perpendicular and parallel to the mean 
magnetic field for magnetic fields of different intensities. 

This work has a broad astrophysical impact. Clusters 
of galaxies is just one of the examples where non-radiative 
heat transfer is essential. This process, however, is impor- 
tant for many regions within galactic interstellar medium, 
e.g. for supernova remnants. 

2. NUMERICAL METHODS 

We use a 3rd-order hybrid essentially non-oscillatory 
(ENO) upwind shock-capturing scheme to solve the ideal 
MHD equations. To reduce spurious oscillations near 
shocks, we combine two ENO schemes. When variables 
are sufficiently smooth, we use the 3rd-order Weighted 
ENO scheme (Jiang & Wu 1999) without characteristic 
mode decomposition. When the opposite is true, we use 
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the 3rd-order Convex ENO seheme (Liu & Osher 1998). 
We use a three-stage Rungc-Kutta method for time inte- 
gration. We solve the ideal MHD equations in a periodic 
box: 

ap/at + v-(pv) = o(i) 

dw/dt + V • Vv + />-^V(aV) - (V x B) x B/47rp = f (2) 

dB/dt-V X (v X B) = 0(3) 

with V • B = and an isothermal equation of state. Here f 
is a random large-scale driving force, p is density, v is the 
velocity, and B is magnetic field. The rms velocity Vturb 
is maintained to be approximately unity, so that v can 
be viewed as the velocity measured in units of the r.m.s. 
velocity of the system and B/y'47r/9 as the Alfven velocity 
in the same units. The time t is roughly in units of the 
large eddy turnover time (^ hnj / Vturh) and the length in 
units of linj , the scale of the energy injection. The mag- 
netic field consists of a uniform background field and a 
fluctuating field: B = Bq -|- b. 

We use a passive scalar to trace thermal particles. 
We inject a passive scalar with a Gaussian profile: 

ip{x,t = to)(x exp-(^-^«)'/'^o , (4) 

where ao= 1/16 of a side of the numerical box and xq lies 
at the center of the computational box. The value of cto 

ensures that the scalar is injected in the inertial range of 
turbulence. The energy injection scale {linj) is ~ 1/2.5 of 
a side of the numerical box. The scalar field follows the 
continuity equation 

di^/dt + V ■ (ipv) = 0. (5) 

We are mainly concerned with time evolution of ctj (i=x, 
y, and z): 



<T? = 



(6) 



where Xi = [/ xtp(x,t)(Px/ J tp(x,t)(Px'\. Common wis- 
dom was that the mean magnetic field suppresses diffusion 
in the direction perpendicular to it. If this is the case, we 
expect to see a± < (T|| . Otherwise, we will get aj_ ^ f7|| . 

We inject passive scalars after turbulence is fully 
developed. Fig. 1(a) shows when we inject the 
passive scalars. For the hydrodynamic run with 
Ms (sonic Mach number) = 0.3 and 192^ grid points 
(thick solid line), we inject passive scalars 5 times. The 
injection times are marked by arrows. We also mark 
the injection times by arrows for the MHD run with 
Va{= Bq/^JWp) = 1, Ms = 0.3, and 192^ grid points 
(thin solid line iov <V^ > and dashed line for < 6^ >). 

3. THEORETICAL CONSIDERATIONS 

Consider two massless particles in the inertial range. Let 
the separation be I. The separation follows 

dl"^ {I + v,dt)- - [I - v,dtf 

'dt Jt 
where we ignore constants of order unity. Using e ~ ''ii /l-, 
we get 



(7) 



^.Z(eOV3, 



(8) 



where e is the energy injection rate. This leads to 

f/^^lT^{Cnf'\''\t-t,l (9) 



where Cr is Richardson constant. When I ^ Iq, we can 
write 

f = CRe{t-t^f, (10) 

which was first discovered by Richardson (1926). 

Recent direct numerical simulations suggest that Cr ~ 
1. Boffetta & Sokolov (2002) obtained Cr ~ 0.55. Ishi- 
hara & Kaneda (2001) obtained Cr ~ 0.7. 

When we inject a passive scalar field as in equation (4), 
we may write 

a'l^-al" = {C,Y/\.'l\t^t,), 

= C2Vturb{t - to), (11) 

where a = (a^ + <^y + o'z)^^^ and the dimensionless con- 
stant Ci is not necessarily the same as Cr. The constant 
C2 cx {Ci/linj)^^^ has dimension. In this paper, we do not 
attempt to obtain Ci or Cr. Instead, we investigate how 
C2 behaves when wc vary Bq. 

Usually it was considered that MHD turbulence is differ- 
ent from its hydrodynamic counterpart. However, recently 
Cho, Lazarian, & Vishniac (2002) showed that motions 
perpendicular to the local mean fields are hydrodynamic 
to high order. This means that many turbulent processes 
are as efficient as hydrodynamics ones. For example, Cho 
et al. (2002) numerically showed that cascade timescale in 
MHD turbulence follows hydrodynamic scaling relations 
(see also Maron & Goldreich 2001). The similarity be- 
tween magnetized and unmagnetized turbulent flows moti- 
vates us to speculate that turbulent mixing is also efficient 
in MHD turbulence. This is why we may use equation 
(11), which is derived from hydrodynamic turbulence. It 
is worth noting that these facts are consistent with a recent 
model of fast magnetic reconnection in turbulent medium 
(Lazarian & Vishniac 1999). 

4. RESULTS 

In Figure 1(b) and (c), we compare the time evolution 
of (J in hydrodynamic case and in MHD case. In the MHD 
case, the Alfven velocity of the mean field (Va = 1) is 
slightly larger than the rms fluid velocity {vturb ^ 0.7). 
This is so-called subAlfvenic regime. Since Va ~ fturb, 
the turbulence is strong. The results show that turbulent 
diffusion is faster in hydrodynamic case. However, Figure 
1(d) implies that this is due to reduction in velocity. Note 
that Vturb ^ 1 in the hydrodynamic case and Vturb ^ 0.7 
in the MHD case (see Figure 1(a)). 

Figure 1(d) shows that there are good relations between 
fT^/3 and {t — tfj). The slopes correspond to the constant 
C2 in equation (11). The slopes are not very sensitive to 
Va or Ms. 

Figure 1(e) and (f) shows that diffusion rate does not 
strongly depend on the direction of the mean field. 
The validity of equation (7) enables us to write 

^dynamic ^dyn^'inj^turb j (-^^) 

where Cdyn is a constant of order unity. This is the ef- 
fective diffusion by turbulent motions suitable for scales 
larger than Imj. The value of Cdyn remains almost con- 
stant for -Bo's of up to Bq ^ 6B ^ vmrb- The exact value 
of Cdyn is uncertain. In hydrodynamic cases, Cdyn is of or- 
der of ~ 0.3 (see Lesieur 1990 chapter VIII and references 
therein) . 
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5. ASTROPHYSICAL IMPLICATIONS 

We have shown that turbulence motions provide efficient 
mixing in MHD turbulence. In this section, we show that 
this process is as efficient as that proposed by Narayan & 
Medvedev (2001) for some clusters. 

We summarize models of thermal diffusion in Fig. 2. 
In the classical picture, thermal diffusion is highly sup- 
pressed in the direction perpendicular to Bq. Transport 
of heat along wondering magnetic field lines (Narayan & 
Medvedev 2001) partially alleviates the problem. But 
the applicability of Narayan & Mcdvcdcv's model is a 
bit restricted - their model requires strong (i.e. Va = 
Bo/^/A-Kp ^ Vturb) niean magnetic field. In the Galaxy, 
there are strong mean magnetic fields. But, in the ICM, 
this is unlikely. When the mean field is weak, the scales 
smaller than the characteristic magnetic field scale (= Ib) 
may follow the Goldreich & Sridhar model (1995). How- 
ever, this requires further studies. Our turbulent mixing 
model gives the same Kdynamic regardless of magnetic field 
geometry. 

ICM — As we mentioned earlier, Ksp ~ 
1.3 X 1030(A;T/A;ey)^/2(n/10-3cm-3)-icm2sec-i and 

Kdynarmc ~ 3.1x10^° {vturb / 500km / s){k„j /20kpc)cm'^ SBC'''' 

The ratio of the two is of order unity for the ICM: 

I^ICM — ^dynamic /i^sp ~ 0(1). (13) 

To be specific, for Hydra A, /x 0.5 for the inner re- 
gion (R ~ lOOkpc) and /x ~ 8.6 for the very inner region 
(R - lOkpc). For 3C 295, /i - 0.34 for the inner region 
and /i ~ 24 for the very inner region. 

Local Bubble and SNRs — The Local Bubble is a hot 
(T ~ lO'^K; kT ~ 100 eV), tenuous (n - O.OOS/cm^) cav- 
ity immersed in the interstellar medium (Berghofer et al 



1998; Smith & Cox 2001). Turbulence parameters are un- 
certain. We take typical interstellar medium values: linj ^ 
10 pc and vturb ^ 5 km/sec. For these parameters, the ra- 
tio of Kdynamic tO KQp is 

Mm = l^dynamic/ I^Sp ^ 0.05, (14) 

for the inside of the Local Bubble. For the mixing layers, 
it is 

l^mix — l^dynamic / ^Sp ^ 100, (-^^) 

where we take f ~ \/T^ ~ 10^ K, n - ^/n^ ~ O.l/cm^ 
(Begelman & Fabian 1990), - 10*K, n, - l/cm^, 
Th ~ lO^if, and nu ~ 0.008/cm^. We expect similar re- 
sults for supernova remnants since parameters are similar. 

6. CONCLUSION 

We have shown that magnetic fields (either random or 
mean magnetic field of up to equipartition value) do not 
suppress turbulent diffusion processes, which implies that 
turbulent diffusion coefficient has the form Kdyn ~ linjVturb 
in MHD turbulence, as well as in hydrodynamic cases. 
This result has two important astrophysical implications. 
First, in the ICM, this turbulent diffusion coefficient is of 
the same order of the classical Spitzer value. Second, in the 
face of hot and cold media in the ISM (e.g. the boundary 
between the Local Bubble and surrounding warm media), 
this turbulent diffusion coefficient is much larger than the 
classical Spitzer value. 
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(d) (e) (f) 

Fig. 1. — Time evolution of energy density and u. (a) Passive scalars are injected after turbulence is fully developed. Injection times are 
marked by arrows. Thick solid line: < V"^ > in the hydrodynamic run. Thin solid line: < V"^ > in the MHD run. Dashed line: < fe^ > in 
the MHD run. (b) Hydrodynamic run. 192^ grid points. (Mach no.) ~ 0.3. (c) MHD run. 192^ grid points. ~ 0.3. Va = 1. (d) 

(p^l^ — ctq^^) /vturb vs. time, a is normalized by the box-size. Y-values are shifted by 0.3 units for convenience, (e) at (i=x, y, and z) vs. 
time. MHD run with 192^ grid points, Ms ~ 0.3, and Va = 1. (f) MHD run with 216^ grid points, Ms ~ 2.3, and Va = 1. Solid lines = x = 
parallel to Bq; dashed lines = y; dotted lines = z. 




Fig. 2. — Models of thermal diffusion, (a) Classical picture. k± <g; Kgp. (b) Narayan & Medvedev (2001). Wandering of field lines provides 
efficient diffusion {k± ~ i^Sp/^) ™ the direction perpendicular to Bq. But, the model assumes Bq of ~ equipartition value, (c) Turbulent 
diffusion model. Thermal electrons are mixed by turbulent motions, which leads to turbulent diffusion coefficient of K^ynamic. ~ '''turbhnj- 
In many astrophysical situations, this coefficient is comparable with the Spitzer value. The figure is the snapshot of the passive scalar field 
at t~3 from the MHD run described in Figure 1(c); 192^ grid points, Ms ~ 0.3, and Va = I- In the case shown here, the mean field is 
strong and parallel to to the dashed line. In general, mean magnetic fields, weak or moderately strong, do not strongly suppress turbulent 
motions / diffusion. 



